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An  analytical  equation  is  derived  for  total  triple  phase  boundary  length  per  unit  volume  (Ltpb)  in  an 
isotropic  uniform  random  microstructure  of  LSM/YSZ  composite  cathode.  The  equation  is  applicable  to 
YSZ  and  LSM  particles  of  any  convex  shapes  and  size  distributions.  The  equation  explicitly  relates  LTpb  to 
the  shapes,  mean  sizes,  coefficient  of  variation  (a  measure  of  the  spread  in  a  size  distribution)  and  skew¬ 
ness  of  YSZ  and  LSM  particle  populations,  and  volume  fractions  of  YSZ,  LSM,  and  porosity.  The  equation 
is  verified  using  available  experimental  data,  and  compared  with  the  results  of  earlier  simulations  and 
models.  The  parametric  analysis  reveals  that  (1)  non-equiaxed  plate-like,  flake-like,  and  needle-like  YSZ 
and  LSM  particle  shapes  can  yield  substantially  higher  LTpb;  (2)  mono-sized  YSZ  and  LSM  powders  lead 
to  higher  LTpb  as  compared  to  the  powders  having  size  distributions  with  large  coefficient  of  variation; 
(3)  Ltpb  is  inversely  proportional  to  the  mean  sizes  of  YSZ  and  LSM  particles;  (4)  high  value  of  LTpb  is 
obtained  at  the  lowest  porosity  volume  fraction  that  permits  sufficient  connectivity  of  the  pores  for  gas 
permeability;  and  (5)  LTpb  is  not  sensitive  to  the  relative  proportion  of  YSZ  and  LSM  phases  in  the  regime 
of  interest  in  composite  cathode  applications. 
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1.  Introduction  and  background 

The  need  for  clean  and  sustainable  energy  has  stimulated  great 
interest  in  fuel  cells.  The  advantages  of  solid  oxide  fuel  cells  (SOFCs) 
over  other  types  of  fuel  cells  include  high  energy  efficiency  and 
excellent  fuel  flexibility  [1,2].  A  wide  variety  of  fuels  can  be  directly 
utilized  in  SOFCs,  including  hydrocarbon  fuels,  coal  gas,  and  renew¬ 
able  fuels.  To  be  economically  competitive  and  commercially  viable, 
however,  the  performance  of  SOFCs  must  be  further  improved  and 
the  cost  must  be  reduced.  The  performance  of  SOFCs  is  often  limited 
by  the  interfacial  losses  arising  from  the  resistances  to  charge  and 
mass  transfer  along  surfaces;  across  interfaces;  and  through  the 
porous  electrodes,  especially  through  the  cathode  where  oxygen 
reduction  takes  place  [3].  Lanthanum  doped  strontium  mangan- 
ite  (LSM)  is  the  most  widely  used  cathode  for  SOFC  based  on 
yttria-stabilized  zirconia  (YSZ)  electrolyte  because  of  its  excellent 
chemical  and  thermal  compatibility  with  YSZ.  Nonetheless,  the 
catalytic  activity  of  LSM  is  severely  limited  by  its  poor  ionic  conduc¬ 
tivity,  particularly  at  low  operating  temperatures.  Consequently,  in 
the  LSM/YSZ  cathodes,  the  electrochemically  active  sites  for  oxy¬ 
gen  reduction  reactions  are  mostly  at  the  triple  phase  boundaries 
(TPB),  which  are  the  lineal  regions  common  to  LSM,  YSZ,  and  air. 


*  Corresponding  author.  Tel.:  +1  404  894  2887;  fax:  +1  404  894  9140. 
E-mail  address:  arun.gokhale@mse.gatech.edu  (A.M.  Gokhale). 

0378-7753 /$  -  see  front  matter  ©  2009  Elsevier  B.V.  All  rights  reserved. 
doi:10.1016/j.jpowsour.2009.05.012 


Therefore,  an  important  strategy  to  increase  the  catalytic  activity  in 
the  cathodes  at  low  temperatures  is  to  increase  the  number  of  such 
active  sites  via  an  increase  in  the  total  length  of  the  TPB  per  unit 
volume  of  microstructure  (Ltpb)  through  rational  design  of  porous 
composites  containing  LSM  and  YSZ  [4-7].  The  performance  of  such 
composite  cathodes  is  microstructure  sensitive  [8-12]. 

The  LSM/YSZ  composite  cathodes  are  typically  fabricated  using 
powder  processing  techniques  involving  printing  or  spraying  and 
sintering  of  a  powder  mix  consisting  of  YSZ  and  LSM  powders. 
It  has  been  reported  that  the  microstructure  of  the  composite 
cathodes  is  sensitive  to  the  processing  parameters  such  as  the 
mean  sizes  of  the  initial  powders  [10,11]  and  sintering  temper¬ 
ature  [9,12].  The  microstructure  in  turn  affects  the  performance 
of  the  cathode.  One  of  the  key  microstructural  parameters  that 
affects  the  electro-chemical  response  of  porous  LSM/YSZ  cathodes 
is  the  total  length  of  the  LSM-YSZ-porosity  triple  phase  bound¬ 
aries  in  the  three-dimensional  (3D)  microstructure  per  unit  volume 
(i.e.,  Ltpb)  because  02  reduction  and  incorporation  of  O2-  into  the 
electrolyte  take  place  at  these  sites  [8].  Accordingly,  development 
of  quantitative  relationships  among  the  processing  parameters, 
microstructural  geometry,  and  electro-chemical  response  of  com¬ 
posite  cathode  materials  is  vital  to  the  effective  optimization  of 
electrode  performance.  Therefore,  it  is  of  interest  to  develop  the 
quantitative  relationships  that  describe  the  effects  of  the  process 
variables  such  as  morphology,  mean  sizes  and  size  distributions  of 
the  powders  of  YSZ  and  LSM  prior  to  sintering  on  LTpb- 
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The  relationships  between  LTPB  and  the  geometric  characteristic 
of  the  powder  mix  prior  to  sintering  have  been  studied  via  geomet¬ 
ric  computer  simulations  of  the  composite  cathode  microstructure 
[13-21  ].  Such  simulated  microstructures  are  also  useful  for  imple¬ 
mentation  in  the  computational  simulations  of  the  electrochemical 
behavior  and  performance  of  the  composite  cathode  [13,14,16-21  ]. 
The  microstructure  simulations  must  allow  overlaps  among  the  YSZ 
and  LSM  particles  to  create  the  TPB  because  simulations  based  on 
non-overlapping  particles  (so  called  hard-core  simulations)  con¬ 
tain  only  point  contacts  among  the  particles.  Consequently,  the 
simulated  LTPB  depends  on  the  extent  of  the  overlaps  permit¬ 
ted.  Geometric  simulations  of  ITPB  reported  in  the  literature  can 
be  divided  into  two  groups:  (a)  the  simulations  where  the  YSZ 
and  LSM  particles  are  placed  on  a  periodic  lattice  but  the  lattice 
points  are  occupied  by  either  YSZ  or  LSM  in  a  random  fashion  with 
the  probability  proportional  to  their  respective  number  densities 
[13,14,16,20],  and  (b)  the  simulations  involving  random  packing  of 
impenetrable  mono-sized  spheres  of  YSZ  and  LSM  followed  by  a 
fixed  amount  of  dilation  (increase)  of  sphere  diameters  to  create 
particle  overlaps  and  TPB  [15,17,18,21  ].  The  stochastic  nature  of  the 
spatial  distribution  of  the  particle  centers  in  the  real  microstruc¬ 
ture  of  the  composite  cathode  is  not  captured  in  the  simulations 
based  on  a  periodic  lattice  (i.e.,  group  (a)),  whereas  the  simula¬ 
tions  based  on  random  packing  of  YSZ  and  LSM  particles  account 
for  the  stochastic  nature  of  the  corresponding  real  microstructures. 
Earlier  simulations  of  the  TPB  assume  that  the  YSZ  and  LSM  parti¬ 
cles  are  spherical  and  mono-sized.  On  the  other  hand,  in  a  typical 
powder  mix,  distributions  of  sizes  of  YSZ  and  LSM  particles  usually 
exist.  Thus,  the  simulations  do  not  reveal  the  effects  of  the  particle 
shape  and  the  distribution  attributes  such  as  the  variance  and  the 
skewness  of  the  particle  sizes  on  the  LtPB. 

Analytical  models  have  also  been  proposed  for  LTPB  [17,20]. 
These  models  assume  that  (i)  YSZ  and  LSM  particles  are  spherical 
and  mono-sized,  (ii)  all  TPB  are  of  the  same  length,  and  (iii)  there  is 
no  distribution  of  coordination  numbers  present  in  the  microstruc¬ 
ture.  To  compute  LTPB  using  these  models,  it  is  essential  to  know 
the  mean  coordination  number  for  the  particle  populations,  which 
cannot  be  measured  experimentally,  except  via  reconstruction 
of  the  3D  microstructure  [22],  and  is  difficult  to  compute  from 
theoretical  considerations  without  making  numerous  simplifying 
assumptions.  Further,  a  distribution  of  coordination  numbers 
usually  exists  in  the  real  microstructures  having  a  distribution  of 
particle  sizes,  and  there  is  often  a  strong  correlation  between  the 
particle  size  and  the  coordination  number:  larger  particles  have 
higher  coordination  numbers  [22].  In  this  contribution,  a  stochastic 
geometry  [23-25]  based  equation  is  derived  that  relates  LTPB  in 
the  3D  microstructure  of  the  composite  cathode  to  the  parameters 
such  as  particle  shape/morphology;  mean  sizes,  standard  deviation 
and  skewness  of  YSZ  and  LSM  size  distributions;  relative  amounts 
of  YSZ  and  LSM;  and  volume  fraction  of  the  porosity.  The  equation 
is  based  on  a  well-known  stereological  relationship  for  estimation 
of  microstructural  surface  area  in  3D  microstructures  [26-28].  The 
derivation  of  the  equation  is  presented  in  the  next  section,  and 
that  is  followed  by  critical  analysis  of  the  model,  comparison  with 
experimental  data,  and  parametric  studies  for  microstructural 
engineering  to  optimize  ITPB. 

2.  Theoretical  development 

2 A.  Extended  microstructure 

Consider  a  geometric  microstructure  model  comprising  an 
isotropic  uniform  random  (IUR)  collection  of  convex1  particles 


1  A  particle  shape  is  convex  if  a  line  joining  any  (and  all)  two  arbitrary  points  on 
its  surface  lies  inside  the  particle.  Spheres,  ellipsoids,  polyhedrons  with  flat  faces, 


of  YSZ  and  LSM  in  a  3D  microstructural  space,  where  the  spa¬ 
tial  arrangement  of  the  particle  centers  is  given  by  the  Poisson 
process  of  spatial  statistics  [29],  and  the  particles  are  allowed  to 
freely  intersect  and  overlap  with  one  another  depending  on  their 
spatial  locations.  Assume  that  there  are  no  spatial  correlations 
in  the  locations  of  the  particles.  The  space  not  occupied  by  YSZ 
and/or  LSM  phases  is  the  porosity.  Such  a  geometric  microstruc¬ 
ture  model  is  called  an  “extended”  microstructure;  the  concept  is 
widely  used  in  modeling  phase  transformations  [30-35].  The  clas¬ 
sical  Johnson-Mehl-Avrami  equation  used  in  phase  transformation 
theories  [30-32]  utilizes  the  concept  of  “extended”  microstructure. 
The  concept  of  “extended”  microstructure  was  initially  developed 
for  two-phase  microstructures,  but  recently  it  has  been  modified  to 
model  microstructures  containing  three  or  more  phases  [33-35]. 

Clearly,  in  a  real  microstructure  the  particles  cannot  overlap. 
Therefore,  the  properties  of  an  extended  microstructure  (such 
as  volume  fractions  of  phases)  are  overestimates  of  the  prop¬ 
erties  of  the  corresponding  “real”  microstructure.  Nonetheless, 
the  “real”  microstructure  model  (i.e.,  microstructure  model  for 
composite  cathode)  can  be  recovered  from  the  corresponding 
extended  microstructure  by  subtracting  the  overlapped  regions 
from  it.  Mathematically,  each  convex  particle  in  an  extended 
microstructure  can  be  considered  as  a  convex  set,  and  an  extended 
microstructure  can  be  treated  as  an  IUR  ensemble  of  convex  sets. 
Therefore,  the  corresponding  “real”  microstructure  is  simply  the 
union  of  these  convex  sets,  which  is  amenable  to  an  analytical  treat¬ 
ment  using  Boolean  algebra  of  convex  sets  [36].  The  microstructural 
properties  of  the  corresponding  “real”  microstructure  obtained  by 
subtracting  the  overlapped  regions  in  the  extended  microstructure 
can  be  computed  analytically  from  the  properties  of  the  extended 
microstructure  [36-38].  Let  (0y)ex  and  (0L)ex  be  the  volume  frac¬ 
tions  of  YSZ  and  LSM  in  our  extended  microstructure.  As  the 
particles  of  YSZ  and  LSM  are  allowed  to  overlap  in  the  extended 
microstructure,  (0Y)ex  and  (0L)ex  are  not  equal  to  their  volume  frac¬ 
tions  0Y  and  Op  in  the  corresponding  real  microstructure,  but  they 
are  related  as  follows  [33-35]. 


(1  +a)0L  =  1  -  exp{-(l  +a)(0L)ex} 

(i) 

Ot  +  1  n  .  f  1+Q.„.  "1 

a  0Y  =  1-eXp{  a  (0yM 

(2) 

where, 

(^Y)ex 

(^L)ex 


(3) 


The  parameter  a  is  the  relative  proportion  of  YSZ  and  LSM  phases 
in  the  microstructure.  In  a  microstructure  containing  YSZ,  LSM,  and 
porosity,  the  sum  of  the  real  volume  fractions  of  YSZ,  LSM,  and 
porosity  must  be  equal  to  one2.  Therefore, 


Oy  +  Op  +  Op  =  \ 


(4) 


In  Eq.  (4),  Op  is  the  volume  fraction  of  the  porosity  in  the  real 
microstructure.  Combining  Eqs.  (1  )-(4)  leads  to  the  following  equa¬ 
tion  for  the  product  of  the  extended  volume  fractions  of  YSZ  and 
LSM.  This  key  equation  will  be  needed  subsequently. 

(0y)ex(0L)ex  =  ; — [ln(6»P  )]2  (5) 

(1  +af 


cylinders,  cones,  plates,  needles,  etc.  are  all  convex. 

2  Note  that  in  the  present  work,  the  volume  fraction  of  a  phase  is  equal  to  the  total 
volume  occupied  by  that  phase  in  the  microstructural  space  divided  by  the  total 
volume  of  the  microstrucutral  space  (i.e.,  specimen  volume).  This  definition  is  dif¬ 
ferent  from  the  one  used  in  the  earlier  papers  on  modeling  of  cathode  microstructure 
[15,17,42]. 
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2.2.  Relationship  between  triple  phase  boundary  lengths  in  real 
and  extended  microstructures 


In  the  extended  microstructure,  the  YSZ  and  LSM  particles  are 
permitted  to  overlap  and  intersect.  The  lines  of  intersection  of 
the  YSZ  and  LSM  particles  are  the  lineal  regions  common  to  YSZ, 
LSM,  and  porosity,  and  therefore,  they  are  the  TPB  of  interest  in 
the  extended  microstructure.  Let  (LTPB)ex  be  the  total  length  of 
these  lines  of  intersection  in  the  extended  microstructure  per  unit 
volume,  and  let  LTpb  be  their  total  length  in  the  corresponding 
real  microstructure  per  unit  volume.  In  general,  LTPB  is  not  equal 
to  (LTPB)ex  because  not  all  segments  of  the  TPB  in  the  extended 
microstructure  contribute  to  the  LTPB  in  the  real  microstructure: 
only  those  triple  phase  boundary  line  segments  are  present  in  the 
real  microstructure  that  are  not  located  in  the  space  already  occu¬ 
pied  by  other  YSZ  and/or  LSM  particles,  i.e.,  those  that  are  located 
in  the  space  occupied  by  the  porosity.  In  an  IUR  microstructure,  the 
probability  that  a  randomly  located  infinitesimal  line  element  falls 
in  the  porosity  phase  is  precisely  equal  to  the  volume  fraction  of  the 
porosity  [23-28],  and  the  probability  is  the  same  for  all  such  line 
elements  and  independent  of  the  location.  Therefore,  the  fraction 
of  the  total  length  of  the  TPB  in  the  extended  microstructure  that 
is  present  in  the  corresponding  real  microstructure,  LTPB/(LTPB)ex,  is 
precisely  equal  to  the  volume  fraction  of  the  porosity  phase,  Op,  i.e., 


£tpb 

(f'TPB  )ex 


=  Op 


or, 


(6) 


LTpb  =  (LtPB  )ex^P 


(7) 


It  remains  to  derive  an  expression  for  (LiPB)ex  in  terms  of  the 
geometric  characteristics  of  YSZ  and  LSM  particles  to  complete  the 
derivation. 


2.3.  Relationship  between  LTPB  and  geometric  attributes  of  YSZ, 

LSM,  and  porosity 

Numerous  geometric  attributes  of  3D  microstructures  can  be 
statistically  estimated  from  the  measurements  performed  on  lower 
dimensional  manifolds  such  as  random  two-dimensional  (2D)  sec¬ 
tions  through  the  3D  microstructure  using  classical  stereological 
relationships  [23-28].  For  example,  volume  fractions  of  the  phases, 
total  surface  areas  of  microstructural  surfaces  per  unit  volume,  and 
integral  mean  curvature  of  surfaces  in  a  3D  microstructure  can  be 
estimated  via  unbiased  sampling  of  the  3D  microstructure  using 
planes  or  surfaces  as  sampling  probes  [23-28].  Consider  estimation 
of  total  surface  area  Si  of  surfaces  of  interest  in  the  3D  microstruc¬ 
ture  of  a  specimen  of  volume  E.  Suppose  this  3D  microstructure 
is  sampled  in  an  unbiased3  manner  with  another  set  of  probe  sur¬ 
faces.  This  can  be  done  by  placing  the  probe  surfaces  of  uniform 
random  orientations  at  a  large  number  of  uniform  random  loca¬ 
tions  in  the  microstructural  space  of  interest.  Let  Sprobe  be  the  total 
area  of  the  probe  surfaces.  The  intersections  of  the  probe  surfaces 
with  the  microstructural  surfaces  of  interest  create  lines  of  intersec¬ 
tion.  Let  Ltotai  be  the  total  length  of  these  lines.  Stochastic  geometry 
gives  the  following  general  relationship  [23-28,39]. 

•^1  4  Ltotal  (q\ 

=  ^  c -  W 

Zj  jl  Sprobe 


3  Unbiased  sampling  implies  that  sampling  at  all  locations  with  probes  having  all 
possible  orientations  is  equally  likely.  Therefore,  each  microstructural  feature  has 
the  same  probability  of  being  included  in  the  statistical  sample. 


or, 

Si  _  4  Ltotal/ £  .  . 

^  Sprobe/^ 

Eqs.  (8)  and  (9)  are  applicable  to  any  3D  microstructure  and  the 
probe  surfaces  of  any  geometry  (for  example,  probes  can  be  sur¬ 
faces  of  ellipsoids,  polyhedrons,  or  planes).  The  only  requirement 
is  that  the  sampling  must  be  unbiased,  which  implies  that  the  sur¬ 
faces  of  the  probes  and  the  surfaces  of  interest  must  have  uniform 
random  orientations  and  locations  with  respect  to  one  another.  This 
requirement  is  satisfied  in  our  extended  microstructure  having  uni¬ 
form  random  orientations  and  locations  of  YSZ  and  LSM  particles. 
Now,  consider  estimation  of  the  total  length  of  the  TPB  per  unit  vol¬ 
ume  (LTPB)ex  in  the  extended  microstructure  using  Eq.  (9).  Consider 
a  thought  experiment  where  the  total  surface  area  of  YSZ  particle 
surfaces  per  unit  volume  (Syz)ex  is  to  be  estimated  (i.e.,  Si/27  in  Eq. 
(9))  using  LSM  particle  surfaces  as  probes  via  application  of  Eq.  (9). 
Let  (SuvOex  be  the  total  area  of  the  LSM  surfaces  per  unit  volume 
(i.e.,  probe  surface  area  per  unit  volume,  SprobeID).  Intersections  of 
YSZ  and  probe  LSM  surfaces  are  the  TPB  whose  total  length  per  unit 
volume  is  (LTPB)ex  (i.e.,  Ltotai/X’).  For  isotropic  uniform  random  YSZ 
and  LSM  particles,  Eq.  (9)  gives  the  following  result4. 


(LrPB)ex 

(SLM)ex  - 
or, 


(Svz)ex  =  ~ 


(10) 


(LTPB)ex  =  ^(5yz)ex(^LM)ex 


(11) 


Combining  Eqs.  (7)  and  (11)  gives  the  following  expression  for 
the  triple  phase  boundary  length  per  unit  volume  LTPB  in  the  corre¬ 
sponding  real  microstructure. 


LtPB  =  ^  ^p(^YZ  )ex(^LM  )ex 


(12) 


Let  (Sy)  and  (Sl)  be  the  mean  values  of  the  surface  areas  and  let 
(Vy)  and  (Vl)  be  the  mean  volumes  of  the  YSZ  and  LSM  particles, 
respectively.  Let  NY  and  NL  be  the  number  of  YSZ  and  LSM  particles 
per  unit  volume,  respectively.  Therefore, 


(Syz)ex  =  (Sy)NY  (13) 

(SuvOex  =  (Sl)Nl  (I4) 

(0Y)ex  =  (W>NY  (15) 

(^L)ex  =  (^L)NL  (16) 


Combining  Eqs.  (5)  and  (13)-(16)  yields  the  following  result. 

(SYz)ex(5 *LM)ex  =  lyilyL  ^  2  lln(#p)]2  07) 

(yy)(Yl)  (1  +a) 

Combining  Eqs.  (12)  and  (17)  gives  the  following  result  for  (LTPB ). 


LtPb  =  ?- — — “2  {$p[ln(0P)]2} 
4  (1  +  a)2  1  J 


'  (Sy)(Sl)  ' 
-<Vy)<Vl>_ 


(18) 


Eq.  (18)  is  applicable  to  YSZ  and  LSM  particles  of  any  convex 
shapes  and  any  size-shape  distributions,  provided  that  the  angular 
orientations  and  locations  of  the  YSZ  and  LSM  particles  are  uniform 
random.  The  YSZ  and  LSM  particles  need  not  be  of  the  same  con¬ 
vex  morphology  or  have  the  same  size-shape  distribution;  the  only 


4  If  the  sampling  surfaces  are  planes,  Eq.  (10)  reduces  to  the  following  stereological 

relationship  [23-28,39],  which  is  used  to  estimate  the  total  microstructural  surfaces 

of  interest  per  unit  volume  S-i/Z  from  the  measurements  of  the  total  length  of  the 

boundary  traces  in  the  random  2D  sections  per  unit  area  IA  [23-25,39]. 
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requirement  is  that  each  particle  must  be  convex.  For  example,  YSZ 
particles  can  be  spheres  and  LSM  particles  can  be  plate  shaped.  The 
result  is  valid  for  any  porosity  volume  fraction  0P,  and  for  any  value 
of  the  ratio  of  the  amounts  of  YSZ  and  LSM,  a.  It  is  convenient  to 
write  Eq.  (18)  in  the  following  form. 


i-TPB  =  jFl(a)F2(0p) 


'  (Sy)(Sl)  ' 
.(WXVl). 


(19) 


where 


F,(a)  = 


a 

(1  +a)2 


and, 


(20) 


F2(0p)  =  0P[ln(0P)]2  (21) 

In  Eq.  (19),  the  term  F2{6P)  explicitly  brings  out  the  dependence 
of  Ltpb  on  the  porosity  volume  fraction  0Pt  whereas  the  term  F^(a) 
captures  the  dependence  of  LTPB  on  the  relative  proportions  of  YSZ 
and  LSM,  a ,  explicitly.  Recall  that  a  is  equal  to  OyIOi  (see  Eq.  (3)). 
Therefore,  the  factors  F^(a)  and  F2{0P)  completely  determine  how 
Ltpb  varies  with  the  volume  fractions  of  the  three  phases,  YSZ,  LSM, 
and  porosity.  The  dependence  of  LTPB  on  the  morphology,  the  mean 
sizes,  and  size-shape  distributions  of  the  YSZ  and  LSM  particles  is 
contained  in  the  last  term  in  the  square  bracket  in  Eq.  (20),  i.e., 
(Sy)(Sl)/(Vy)(Vi).  The  mean  volumes  and  surface  areas  of  YSZ  and 
LSM  particles  can  be  expressed  as  follows. 


(Sy)=K,y(D2)  (22) 

<SL>=K1L(D2>  (23) 

(Vy)=K2y(D3)  (24) 

(Vl)=K21(DI)  (25) 


In  the  above  equations,  I< 1Y,  /<tl,  /<2y,  f^2L  are  the  shape  factors 
that  depend  on  the  shape(s)  of  the  YSZ  and  LSM  particles.  (D2)  is  the 
mean  value  of  the  square  of  the  YSZ  particle  sizes  (i.e.,  the  second 
moment  of  the  size  distribution),  and  (D2)  is  the  mean  value  of 
square  of  the  LSM  particle  sizes5.  Similarly,  (D5)  is  the  mean  value 
of  the  cube  of  the  YSZ  particle  sizes  (i.e.,  the  third  moment  of  the 
size  distribution),  and  (D5)  is  the  mean  value  of  cube  of  the  LSM 
particle  sizes.  Combining  Eq.  (19)  with  the  Eqs.  (22)-(25)  gives  the 
following  result. 


Itpb  =  jF1(a)F2(0p)F3(K)^|^  (26) 

where 

K2Yk2L 

F3(I< )  depends  only  on  the  shapes  of  the  YSZ  and  LSM  particles. 
For  spherical  YSZ  and  LSM  particles,  F3(I< )  is  equal  to  36  and  the 
diameters  are  the  size  parameters.  It  remains  to  deconvolute  the 
effects  of  the  mean  particle  sizes,  spreads,  and  skewness  of  the  size 
distributions  of  YSZ  and  LSM  on  the  LTpb  that  are  contained  in  the 
term  (D2)  (D2) /  (D5)  (D5)  in  Eq.  (26).  Let  D  represent  the  particle  sizes 
in  a  size  distribution.  The  arithmetic  mean  particle  size  (D),  the 
variance  cr2,  the  coefficient  of  variation  CV,  and  the  skewness  y  of 
a  size  distribution  function/(D)  are  defined  as  follows  [40,41  ]. 


(D) 


Df(D)dD 


(28) 


5  In  general,  the  mean  of  the  square  of  the  particles  sizes  is  not  equal  to  the  square 
of  the  mean  value  of  the  particle  sizes,  except  when  all  particles  are  of  the  same  size. 
Similarly,  the  mean  of  the  cube  of  the  particles  sizes  is  not  equal  to  the  cube  of  the 
mean  value  of  the  particle  sizes,  except  when  all  particles  are  of  the  same  size. 


a2=  [D  -  (D)]  f(D)dD  =  (D2)  -  (D)2 


(29) 


cv=  2- 

(D) 


Y=~^ 

(7 


hi 


[D  -  (D)]3f(D)dD 


Combining  Eqs.  (29)-(31)  gives  the  following  results. 

<D2>  =  (D)2(l  +CV2) 

(D3)  =  <D)3(1  +  3CV2  +  yCV3) 


(30) 

(31) 

(32) 

(33) 


Substituting  Eqs.  (32)  and  (33)  into  Eq.  (26)  leads  to  the  following 
result. 


1-tpb  =  |F1(a)F2(0P)F3(lC)F4(CV,  y) 


1 

(Dy)(Dl) 


where 


f4(CV,  y)  = 


_ (1+CV2)(1+CV2) _ 

(1  +  3CV2  +  KyCV?)(1  +  3CV2  +  yLCV^) 


(34) 


(35) 


In  Eq.  (34),  CVY  and  CVL  are  coefficient  of  variation,  and  yY  and 
y l  are  skewness  parameters  of  the  YSZ  and  LSM  powder  size  dis¬ 
tributions,  respectively.  Note  that  F4(CV,y)  depends  only  on  the 
coefficient  of  variation  and  the  skewness  of  the  size  distributions, 
and  it  explicitly  brings  out  the  effect  of  these  parameters  on  LTPB. 

Eq.  (34)  predicts  that  for  given  volume  fractions  of  YSZ,  LSM, 
and  porosity,  and  the  coefficient  of  variation  and  the  skewness  of 
the  YSZ  and  LSM  size  distributions,  and  particle  morphologies,  LTpb 
is  inversely  proportional  to  the  mean  size  of  YSZ  particles,  (DY), 
and  the  mean  size  of  LSM  particles,  (Dl).  Therefore,  the  finer  the 
mean  particle  sizes  of  YSZ  and/or  LSM,  the  higher  is  the  total  length 
of  the  TPB  per  unit  volume,  LTPB.  Similar  trend  has  been  predicted 
by  earlier  simulations  and  analytical  models  for  LTPB  in  compos¬ 
ite  cathodes  [15,17,42].  Nonetheless,  earlier  analytical  treatments  of 
Ltpb  were  based  on  an  assumption  that  all  YSZ  and  LSM  particles  are 
of  the  same  size,  and  therefore,  did  not  capture  the  separate  effects 
of  different  YSZ  and  LSM  mean  sizes  on  the  LTpb-  The  present  result 
shows  that  the  LTpb  can  be  increased  either  by  decreasing  the  mean 
size  of  YSZ  particles,  or  LSM  particles,  or  by  decreasing  the  mean 
sizes  of  both  the  particle  populations,  and  it  explicitly  captures  the 
dependence  of  LTpb  on  the  mean  particle  sizes.  Eq.  (34)  also  explic¬ 
itly  brings  out  the  effects  of  other  geometric  characteristics  on  LTpb 
through  the  functions  F\  (a),  F2{0P)}  F3{I< ),  and  F4(CV,y),  which  facil¬ 
itates  the  parametric  studies.  Analysis  of  the  results,  comparison 
with  experimental  data,  and  parametric  studies  for  microstructural 
engineering  to  optimize  LTpb  are  presented  in  the  next  section. 


3.  Results  and  discussion 

A  stochastic  geometry  based  approach  is  presented  for  geomet¬ 
ric  modeling  of  total  triple  phase  boundary  length  per  unit  volume 
Ltpb  in  the  microstructure  of  composite  SOFC  cathodes.  The  model 
is  based  on  the  assumptions  that  the  YSZ  and  LSM  particles  are 
convex,  they  have  uniform  random  angular  orientations  and  spa¬ 
tial  locations,  and  there  are  no  spatial  correlations  in  their  locations. 
The  model  permits  the  particles  to  freely  overlap  and  intersect  in  the 
extended  microstructure.  The  approach  enables  tractable  analytical 
treatment  of  the  geometric  problem,  and  leads  to  closed  form  ana¬ 
lytical  solution  that  relates  LTpb  to  numerous  geometric  parameters 
including  volume  fractions  of  the  phases;  morphologies/shapes  of 
particles;  and  the  mean  sizes,  the  coefficient  of  variation,  and  the 
skewness  of  YSZ  and  LSM  size  distributions.  Although  Eq.  (34)  is 
derived  in  the  context  of  the  microstructures  of  SOFC  composite 
cathodes,  it  is  equally  applicable  to  the  microstructures  of  SOFC 
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anodes,  or  to  any  three  phase  uniform  random  isotropic  microstruc¬ 
ture.  Unlike  the  earlier  analytical  models  [17,42],  to  calculate  the 
Ltpb  using  Eq.  (34),  it  is  not  necessary  to  know  the  mean  coor¬ 
dination  number  of  the  particles,  which  is  difficult  to  measure 
experimentally  (because  reconstruction  of  3D  microstructure  is 
required  [22]),  and  is  difficult  to  compute  from  theoretical  consider¬ 
ations  without  making  numerous  simplifying  assumptions.  Further, 
earlier  models  and  simulations  assumed  that  YSZ  and  LSM  particles 
are  mono-sized  spheres,  and  therefore,  did  not  reveal  the  effects 
of  the  morphology/shape  of  the  YSZ  and  LSM  particles  and  the 
powder  characteristics  such  as  the  coefficient  of  variation  (a  mea¬ 
sure  of  spread  in  the  particle  sizes)  and  the  skewness  of  the  size 
distributions  on  LTPB. 

Eq.  (34)  is  developed  for  prediction  of  LTPB  in  powder-processed 
composite  cathodes  fabricated  via  sintering  of  a  powder  mix  con¬ 
taining  YSZ  and  LSM  powder  particles.  It  may  have  to  be  modified 
for  the  composite  cathodes  fabricated  using  other  processing 
routes.  It  is  important  to  point  out  that  the  present  model  pre¬ 
dicts  the  total  (active  plus  inactive)  triple  phase  boundary  length. 
Some  TPB  are  not  topologically  connected  to  the  electrical  con¬ 
ducting  paths  and  gas  permeating  pores  in  the  microstructure,  and 
therefore,  do  not  participate  in  the  electrochemical  processes.  In  a 
subsequent  contribution  it  will  be  shown  that  beyond  the  percola¬ 
tion  thresholds,  about  90%  of  the  total  triple  phase  boundary  length 
is  connected  to  the  electrical  conducting  paths  and  gas  permeating 
pores  in  the  microstructure  [51  ]. 

All  TPB  that  are  connected  to  the  electrical  conductive  paths  and 
gas  permeating  pores  also  may  not  be  electrochemically  active.  The 
electrochemical  activity  at  the  TPB  depend  on  the  availability  of 
the  electro-active  species  involved  in  the  electrochemical  reactions 
(e.g.,  V0 ••,  h*,  and  02)  at  or  near  the  TPB,  which  is  dictated  by  V## 
transport  through  YSZ,  h ••  transport  through  LSM,  and  02  transport 
through  the  pores  of  the  electrode.  For  a  typical  YSZ/LSM  cathode 
with  sufficient  porosity  for  02  transport,  the  electrochemical  activ¬ 
ity  induced  by  a  triple  phase  boundary  often  diminishes  with  the 
distance  from  the  interface  between  the  YSZ/LSM  composite  cath¬ 
ode  and  the  YSZ  electrolyte  because  the  ionic  conductivity  of  YSZ 
is  orders  of  magnitudes  smaller  than  the  electronic  conductivity  of 
LSM.  Therefore,  only  the  triple  phase  boundary  segments  within  a 
certain  effective  membrane  thickness  (which  is  less  than  the  geo¬ 
metric  thickness  of  the  composite  cathode)  that  are  connected  to 
the  electrical  conductive  paths  and  gas  permeating  pores  are  elec¬ 
trochemically  active  [17,50].  Nonetheless,  it  can  be  said  that  in  a 
composite  cathode  having  isotropic  uniform  random  microstruc¬ 
ture,  an  increase  in  the  LTPB  is  expected  to  lead  to  an  increase  in  the 
electrochemically  active  triple  phase  boundary  length  per  unit  vol¬ 
ume.  Therefore,  the  microstructural  engineering  to  optimize  the 
LtPB  using  the  present  approach  should  also  lead  to  an  increase 
in  the  electrochemically  active  triple  phase  boundary  length  per 
unit  volume,  which  is  likely  to  improve  the  performance  of  the 
cathode. 

The  present  model  is  applicable  only  to  the  microstructural  vol¬ 
ume  segments  that  are  at  least  of  the  size  of  the  smallest  statistical 
representative  volume  element6  (RVE)  of  the  corresponding  infinite 
global  microstructure.  This  is  true  of  any  theoretical  model  or  simu¬ 
lation.  Therefore,  the  present  results  may  not  be  applicable  to  very 
thin  composite  cathode  membranes  (i.e.,  thickness  on  the  order 
of  the  particle  size).  In  a  subsequent  contribution,  we  will  exam¬ 
ine  at  what  ratio  of  membrane  thickness  to  mean  particle  size  the 
microstructural  segment  becomes  sufficiently  large  to  be  treated  as 
RVE. 


6  A  representative  volume  element  of  a  microstructure  is  such  that  if  different 
volume  elements  of  that  size  are  chosen  at  random  anywhere  in  the  corresponding 
infinite  microstructure  all  of  them  will  have  the  same  properties. 


The  input  size  distribution  data  (mean  sizes,  coefficient  of  vari¬ 
ation,  and  skewness)  needed  for  the  present  model  pertain  to  the 
YSZ  and  LSM  particles  present  in  the  powder  mix  after  ball  milling 
(or  after  any  attrition  process  used  to  obtain  a  homogeneous  pow¬ 
der  mix)  and  prior  to  sintering  (i.e.,  not  size  distributions  in  the 
unmixed  initial  powders  of  YSZ  and  LSM)  because  operations  like 
ball  milling  fragment  the  particles  and  alter  the  size  distributions. 
Further,  if  the  powder  particles  agglomerate  during  any  of  the  pow¬ 
der  processing  steps  then  the  distribution  characteristics  of  the 
agglomerated  powders  must  be  used  in  Eq.  (34). 

The  present  approach  is  purely  geometric:  it  does  not  incorpo¬ 
rate  the  physics  of  the  sintering  process.  Therefore,  in  its  present 
form,  the  model  does  not  account  for  coarsening  of  the  microstruc¬ 
ture  during  sintering  that  can  decrease  the  LTPB.  In  the  case  of 
LSM/YSZ  composite  cathodes,  the  sintering  is  typically  carried  out 
at  temperatures  below  1200  °C  to  minimize  such  coarsening  so  that 
the  fine  microstructure  (and  higher  ITPB)  can  be  retained  after  sin¬ 
tering  [11,12],  and  in  such  cases,  Eq.  (34)  can  predict  the  ITPB  in  a 
precise  manner. 

It  is  important  to  point  out  that  the  earlier  modeling  and  sim¬ 
ulation  efforts  were  also  either  geometric  or  involved  substantial 
simplifying  assumptions  concerning  the  microstructural  geometry 
and  the  physics  of  the  sintering  process.  This  is  because  a  geometri¬ 
cally  general  and  rigorous  sintering  physics  based  modeling  of  LTPB 
is  too  complex  to  be  analytically  tractable  or  computationally  fea¬ 
sible.  Consequently,  it  is  of  interest  to  compare  the  predictions  of 
the  present  model  with  the  experimental  data  on  LTPB  and  earlier 
models  and  simulations. 


3.1  Comparison  with  experimental  data 


To  verify  the  present  model  (Eq.  (34)),  the  following  experimen¬ 
tal  data  are  needed:  (i)  experimentally  estimated  LTPB  either  via 
reconstruction  of  3D  microstructure  or  by  using  the  stereological 
techniques,  (ii)  volume  fractions  of  YSZ,  LSM,  and  porosity,  and  (iii) 
mean  sizes,  coefficient  of  variation,  and  skewness  of  the  YSZ  and 
LSM  particle  size  distributions  after  ball  milling  or  powder  mix¬ 
ing  operations  and  prior  to  sintering.  There  have  been  only  two 
experimental  measurements  of  LjPB  in  composite  SOFC  cathode 
microstructures  reported  in  the  literature  [47,49],  and  one  study 
has  reported  LTPB  data  in  SOFC  anode  microstructure  [48].  In  all  of 
these  cases,  the  particle  size  distribution  data  needed  to  verify  Eq. 
(34)  have  not  been  reported.  Nonetheless,  Eq.  (12)  (from  which  Eq. 
(34)  has  been  obtained  without  any  additional  assumptions)  can 
be  cast  into  a  form  that  can  be  verified  using  the  available  data.  For 
an  isotropic  uniform  random  microstructure,  it  can  be  shown  that 
[43,44]: 


SYP  =  0P(SY)ex  (36) 

and, 

SLP  =  0?(Sl)ex  (37) 


SYP  and  SlP  are  the  total  surface  areas  of  YSZ/porosity  surfaces 
and  LSM/porosity  surfaces  per  unit  volume  of  specimen  in  the  real 
microstructure,  which  can  be  experimentally  measured.  Therefore, 
combining  Eqs.  (12),  (36)  and  (37)  leads  to  the  following  result. 


^tpb  = 


it  SypSlP 
4  "" 0p~ 


(38) 


An  equation  analogous  to  Eq.  (38)  was  reported  earlier  in  the 
context  of  microstrucutral  evolution  during  phase  transformations 
[45].  It  is  important  to  point  out  that  Eq.  (38)  is  applicable  to 
any  isotropic  uniform-random  microstructure:  its  applicability  is 
not  limited  to  the  microstructures  of  powder  processed  composite 
cathodes.  Nonetheless,  the  use  of  Eq.  (38)  requires  the  experimen¬ 
tal  microstrucutral  data,  and  therefore,  unlike  Eq.  (34),  it  cannot 
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predict  Ltpb  in  an  a  priori  manner  from  initial  raw  material  charac¬ 
teristics  (for  example,  size  distributions  of  initial  powder  particles) 
before  fabrication  of  the  composite  cathode.  Wilson  et  al.  [46] 
reconstructed  3D  microstructure  of  SOFC  anode  containing  YSZ, 
Ni,  and  porosity,  and  measured  LTPB,  porosity  volume  fraction  0P , 
and  total  surface  areas  of  YSZ-porosity,  SYp,  and  Ni-Porosity,  SNP 
per  unit  volume.  As  the  anode  contains  Ni  instead  of  LSM,  in  Eq. 
(38),  SLp  has  to  be  replaced  by  SNP.  Wilson  et  al.  reported  ITPB 
equal  to  4.28  fxm  p,m-3,  and  0P,  Syp,  and  Snp  equal  to  0.195,  0.5 
and  1.9  p,m2  pm-3,  respectively.  Substituting  these  data  in  Eq.  (38) 
yields  the  value  of  ITPB  equal  to  3.83  pm  pm-3,  which  is  in  good 
agreement  with  the  experimental  value  of  4.28  pm  pm-3. 

Zhang  et  al.  [49]  estimated  the  ITpB  using  atomic  force 
microscopy  and  stereological  techniques  in  a  composite  LSM/YSZ 
cathode7  to  be  equal  to  10.8  ±  1.15  pm  pm-3.  The  experimentally 
measured  volume  fractions  of  YSZ  and  LSM  were  (0.28  ±0.027) 
and  (0.26  ±0.018),  respectively,  and  the  experimentally  measured 
SYp  and  SLp  were  (3.87  ±0.21)  and  (2.32  ±0.3)  pm2  pm-3,  respec¬ 
tively.  Recall  that  0P  is  equal  to  (1  -0Y-0i).  Substituting  these 
data  into  Eq.  (38)  yields  the  computed  value  of  ITpB  equal  to 
(15.8±4.3)pmpm-3.  The  large  error  bar  in  the  computed  LTPB  is 
because  each  of  the  three  experimentally  measured  parameters  on 
the  right  hand  side  in  Eq.  (38)  has  a  statistical  sampling  error.  Con¬ 
sidering  large  sampling  error  associated  with  computed  LTPB,  there 
is  reasonable  agreement  between  the  computed  and  experimen¬ 
tally  measured  ITPB  values. 

Wilson  et  al.  [47]  performed  3D  microstructure  reconstructions 
of  a  series  of  composite  cathodes  having  different  proportions 
of  YSZ  and  LSM,  and  reported  that  the  ITPB  is  in  the  range  of 
8-10  pm  pm-3  and  it  is  not  sensitive  to  the  proportion  of  YSZ  and 
LSM.  This  trend  can  be  predicted  from  Eq.  (34),  where  the  depen¬ 
dence  of  Ltpb  on  the  relative  proportion  of  YSZ  and  LSM  represented 
by  parameter  a  resides  in  the  function  F\  (a).  When  a  is  in  the  range 
0.5-2. 5,  F^{a)  is  in  the  range  of  (0.225  ±0.0225),  i.e.,  a  variation 
of  ±10%.  Thus,  ITPB  is  not  sensitive  to  relative  proportion  of  YSZ 
and  LSM  when  a  is  in  the  range  0.5-2.5  (the  range  for  practical 
applications),  as  observed  by  Wilson  et  al.  [47]. 

3.2.  Comparison  with  simulations  and  models  reported  in 
literature 


lated  ITpB  obtained  by  Ali  et  al.  [15]  and  the  corresponding  values 
computed  using  Eq.  (34). 

Schneider  et  al.  [17]  have  proposed  the  following  analytical 
equation  for  total  (active  plus  inactive)  LT PB  based  on  random  pack¬ 
ing  of  mono-size  spheres  of  radius  r. 


^tpb  =  \j^~  (~j)  faott  - 

where, 

^  _  0Y  _  a 
* i0  0Y  ± 0i  1+a 


(39) 


(40) 


In  Eq.  (39),  d0  is  the  density  of  the  powder  mix  before  sintering, 
d  is  the  density  after  sintering,  and  Z  is  the  mean  coordination  num¬ 
ber  of  the  particles.  Interestingly,  the  0IO(1  -  0lo)  factor  in  Eq.  (39)  is 
equal  to  a/{  1  +  a)2,  i.e.,  Fi  (a)  in  Eq.  (34),  and  for  mono-sized  YSZ  and 
LSM  having  the  same  radius  r,  Eqs.  (39)  and  (34)  predict  that  LTPB  is 
inversely  proportional  to  the  square  of  the  particle  size.  However, 
Eq.  (34)  is  applicable  to  YSZ  and  LSM  particles  having  any  size  distri¬ 
bution  and  convex  particle  shape(s).  On  the  other  hand,  Schneider 
et  al.  assume  that  YSZ  and  LSM  particles  are  spherical  and  mono¬ 
sized,  and  to  compute  LTPB  using  their  equation  it  is  necessary  to 
assume  a  value  for  the  mean  coordination  number  Zand  the  density 
of  initial  powder  mix  d0.  Further,  Eq.  (39)  leads  to  a  physically  unac¬ 
ceptable  limit  as  d  approaches  1  (i.e.,  fully  dense  material).  In  the 
limit  of  d  ->  1  (i.e.,  porosity  volume  fraction  0P  ->  0)  the  LTPB  must 
approach  zero  because  YSZ-LSM-porosity  TPB  cannot  exist  when 
there  are  no  pores.  However,  Eq.  (39)  predicts  a  non-zero  value  of 
Ltpb  as  d  -*  1,  which  is  physically  untenable.  On  the  other  hand,  the 
LjpB  predicted  by  Eq.  (34)  reaches  the  correct  limit  of  zero  as  the 
volume  fraction  of  porosity  0P  approaches  zero  or  one.  Schneider  et 
al.  computed  the  values  of  LTPB  from  Eq.  (39)  assuming  d0  =  0.5  and 
Z= 6.3.  For  these  values  of  d0  and  Z,  Eq.  (39)  gives  LTPB  in  the  range  of 
1.1 -1.2  pun  pun-3  when  r  =  1  pun  and  D  =  0.75  for  <pio  in  the  range  of 
0.3-0.7.  For  the  same  range  of  0IO,  and  r  and  d  values,  Eq.  (34)  yields 
Ltpb  in  the  range  of  0.7-0.84  pm  pm-3,  which  is  somewhat  lower 
than  that  obtained  from  Eq.  (39).  These  differences  are  probably 
due  to  the  values  of  d0  and  Z  assumed  by  Schneider  et  al.  (some¬ 
what  different  values  of  these  parameters  lead  to  better  agreement) 
and/or  because  Eq.  (39)  is  inaccurate  at  high  values  of  d,  as  it  gives 
physically  unacceptable  nonzero  value  of  LTpB  as  d  approaches  one. 


Ali  et  al.  [15]  performed  simulations  of  composite  cathode 
microstructures  involving  random  packing  of  impenetrable  mono¬ 
sized  spheres  followed  by  10%  dilation  of  the  particle  diameters  to 
create  the  TPB.  The  simulations  predict  that  the  LTpB  decreases  with 
increasing  particle  size,  with  a  trend  similar  to  that  predicted  by 
Eq.  (34).  Nonetheless,  the  simulated  triple  phase  boundary  length 
reported  by  these  researchers  is  the  active  LTpB*  whereas  Eq.  (34) 
predicts  the  total  (active  plus  inactive)  Ltpb-  As  mentioned  earlier, 
beyond  the  percolation  threshold,  about  90%  or  more  of  the  total 
Ltpb  is  active  LTpB  [51],  and  therefore,  there  is  no  significant  differ¬ 
ence  between  the  two.  Consequently,  the  simulated  LTpB  reported 
by  Ali  et  al.  [15]  can  be  compared  with  the  predictions  of  Eq.  (34) 
for  the  porosity  volume  fractions  in  the  range  of  0.3-0.5,  and  rel¬ 
ative  proportions  of  YSZ  and  LSM  (i.e.,  parameter  a  in  Eq.  (34))  in 
the  range  of  0.5-2.5.  Ali  et  al.  [15]  assumed  that  the  YSZ  and  LSM 
particles  are  mono-sized  spheres.  For  these  particle  characteristics, 
in  Eq.  (34),  F3{I<)  =  36  and  F4(CV,y)  =  1.  The  LTpb  computed  from  Eq. 
(34)  using  these  values  of  F3(I< )  and  F4(CV,y)  for  the  porosity  vol¬ 
ume  fractions  of  0.3  and  a  in  the  range  of  0.5-1.5  are  given  in  Table  1 
for  comparison  with  the  simulated  LTpb  values  obtained  by  Ali  et 
al.  [15].  There  is  a  good  agreement  between  the  values  of  simu¬ 


7  The  error  bars  are  statistical  sampling  errors  for  95%  confidence  interval. 


3.3.  Parametric  studies 

Eq.  (34)  predicts  that  for  given  morphologies,  volume  fractions, 
coefficient  of  variation,  and  skewness  of  the  YSZ  and  LSM  particle 
populations,  the  LTpB  is  inversely  proportional  to  the  mean  sizes  of 
YSZ  and  LSM  particles,  (Dy)  and  (DL).  This  trend  is  also  predicted  by 
earlier  simulations  and  analytical  models  [15,17,42].  Nonetheless, 
earlier  analytical  treatments  of  LTpB  assumed  that  YSZ  and  LSM  par¬ 
ticles  are  spherical  and  of  the  same  mono-size,  and  therefore,  did 
not  reveal  the  separate  effects  of  different  YSZ  and  LSM  mean  sizes 
on  the  LTpB.  The  present  result  shows  that  LTpb  can  be  increased  by 
decreasing  the  mean  size  of  YSZ  or  LSM  particles,  or  of  both  YSZ  and 
LSM  particles,  and  it  explicitly  captures  the  dependence  of  LTpb  on 
the  mean  particle  sizes. 

3.3 A.  Effect  of  relative  proportion  of  YSZ  and  LSM  on  LTPB 

The  relative  proportion  of  YSZ  and  LSM  is  specified  by  the  param¬ 
eter  a ,  which  is  equal  to  the  ratio  of  the  volume  fraction  of  YSZ  and 
LSM,  [0y/0L].  In  Eq.  (34),  the  dependence  of  LTpb  on  a  is  expressed 
by  the  function  Fi(a).  In  principle,  a  can  vary  from  zero  to  infin¬ 
ity,  and  at  these  two  limits,  F^a)  and  LTpb,  approach  zero  as  they 
must  for  the  model  to  be  physically  tenable.  Nonetheless,  for  SOFC 
applications,  YSZ,  LSM,  and  porosity  phases  must  topologically  per¬ 
colate.  Therefore,  each  phase  should  have  volume  fraction  in  the 
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Comparison  of  the  TPB  length  reported  in  Ali  et  al.[15]  and  calculated  with  Eq.  (34). 
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range  0.25-0.5  such  that  the  sum  of  the  three  volume  fractions  is 
equal  to  one.  With  these  constraints  on  the  volume  fractions,  the 
parameter  a  is  constrained  in  the  range  0.5-2.5.  Fig.  1  shows  the 
variation  of  F^  ( a )  with  a.  Observe  that  F\  (a)  reaches  the  maximum 
value  of  0.25  when  a  is  equal  to  1,  and  it  varies  only  from  0.225  to 
0.25  as  a  varies  from  0.5  to  2.5.  Consequently,  LTPB  does  not  vary  by 
more  than  10%  in  the  range  of  a  values  of  interest  in  SOFC  compos¬ 
ite  cathode  applications  where  all  three  phases  must  topologically 
percolate. 

3.3.2.  Effect  of  porosity  volume  fraction  on  LTPB 

In  Eq.  (34),  the  dependence  of  LTPB  on  the  porosity  volume 
fraction  0P  is  captured  in  the  function  F2(0P ),  which  is  equal 
to{0P[ln0P]2}.  In  principle,  6P  can  vary  from  zero  to  one.  At  these 
two  limits,  F2{0p)  and  consequently  LTPB,  approach  zero  as  they 
must  for  the  model  to  be  physically  tenable.  Fig.  2  shows  the  vari¬ 
ation  of  F2{0p)  with  0P.  Note  that  F2{0P)  has  a  maximum  value  of 
0.541,  when  0P  is  equal  to  0.135.  Nonetheless,  for  practical  SOFC 
composite  cathode  applications,  YSZ,  LSM,  and  porosity  phases, 
must  percolate.  Therefore,  0P  is  expected  to  be  in  the  range  of 
0.2-0.5.  As  6P  increases  from  0.2  to  0.5,  F2{0P)  decreases  from  0.52 
to  0.24  (i.e.,  by  a  factor  of  about  two).  Thus,  for  SOFC  composite 
cathode  applications,  for  given  values  of  a ,  and  mean  particle  sizes, 
CV,  and  the  skewness  of  the  YSZ  and  LSM  powder  populations,  the 
highest  value  of  ITPB  is  obtained  at  the  lowest  porosity  volume  frac¬ 
tion  that  permits  sufficient  percolation  and  connectivity  of  pores 
for  gas  permeability. 

3.3.3.  Effects  ofCV  and  skewness  on  LTPB 

In  Eq.  (34),  for  given  mean  sizes  of  YSZ  and  LSM  powder  par¬ 
ticles,  the  dependence  of  LTPB  on  the  coefficient  of  variation  CV 
and  the  skewness  of  YSZ  and  LSM  particle  size  distributions  is  con¬ 
tained  in  the  function  F4(CV,y),  given  by  Eq.  (35).  For  mono-sized 
YSZ  and  LSM,  coefficients  of  variation,  CVY  and  04,  and  skew¬ 
ness,  yy  and  y i  are  zero,  and  therefore,  F4(CV,y)  is  equal  to  one. 
It  is  well  known  that  for  any  physically  realizable  size  distribution 
function,  skewness  y>(-2/CV),  and  therefore,  F4(CV,y)  <  1.  Con¬ 
sequently,  for  given  relative  proportion  of  YSZ  and  LSM,  porosity 
volume  fraction,  and  mean  particle  sizes  of  YSZ  and  LSM,  the  LTPB 
has  the  maximum  value  for  mono-size  powders  of  YSZ  and  LSM.  Any 


spread  in  the  size  distribution  (i.e.,  non-zero  CV)  decreases  LtPb-  For 
symmetric  size  distribution  functions  such  as  the  normal  distribu¬ 
tion,  the  skewness  y  is  equal  to  zero.  For  such  populations,  as  CVY 
and  CVL  become  very  large  (strictly  speaking,  as  they  go  to  infinity), 
F4(CV,y)  approaches  1  /9.  Therefore,  LtPb  can  decrease  almost  by  an 
order  of  magnitude  as  coefficients  of  variation  of  YSZ  and  LSM  pow¬ 
der  populations  increase  from  zero  (corresponding  to  mono-size) 
to  infinity.  Thus,  there  is  a  strong  dependence  of  F4(CV,y)  on  CV.  For 
example,  relatively  small  values  CVY  =  CVL  =  0.5  in  normal  distribu¬ 
tions  of  YSZ  and  LSM  powders  reduce  F4  (CV,y)  to  0.51  from  the 
value  of  1.0  for  mono-size  powders  (i.e.,  CVY  =  04  =  0),  which  leads 
to  a  decrease  in  ITPB  approximately  by  a  factor  of  2  as  compared 
to  that  for  the  mono-sized  powders.  Therefore,  mono-size  particle 
populations  (although  YSZ  and  LSM  can  have  different  mono-sizes) 
optimize  LtPb  for  given  mean  sizes,  shapes,  and  volume  fractions  of 
YSZ  and  LSM. 

Skewness  of  a  size  distribution  (see  Eq.  (31))  can  be  negative, 
positive,  or  zero.  Inspection  of  Eq.  (35)  reveals  that  for  given  values 
of  CV  of  YSZ  and  LSM  particles,  F4(CV,y)  is  higher  (and  therefore, 
Ltpb  is  higher)  for  YSZ  and  LSM  populations  that  have  a  negative 
skewness ,  than  for  the  size  distributions  that  have  positive  skewness. 
In  a  distribution  having  negative  skewness,  the  mean  is  lower  than 
the  median  (which  is  lower  than  the  mode).  Therefore,  negative 
skewness  leads  to  a  long  left  tail  in  a  size  distribution.  As  particle 
sizes  cannot  be  negative,  clearly  there  is  a  limit  on  the  extent  to 
which  the  skewness  can  be  negative  in  a  physically  realizable  size 
distribution  function.  It  is  not  clear  if  the  current  powder  processing 
technologies  can  yield  YSZ  and  LSM  powder  populations  with  a 
specified  CV  and/or  skewness  values,  but  if  that  is  possible,  Eqs.  (34) 
and  (35)  provide  a  framework  to  compute  how  much  improvement 
in  the  LTPb  is  possible  via  such  process  designs. 

3.3.4.  Effects  of  powder  particle  shape/morphology  on  LTPB 

Earlier  investigations  on  modeling  and  simulations  of  LTPB 
assumed  that  YSZ  and  LSM  particles  are  spherical,  and  therefore, 
did  not  reveal  the  effects  of  particle  shape(s)  on  ITPB.  In  Eq.  (34), 
the  effect  of  particle  shape(s)  on  ITPB  is  contained  in  the  function 
F3(/C),  defined  in  Eq.  (27).  For  spherical  YSZ  and  LSM  particles,  F3{I<) 
is  equal  to  36.  For  given  size  distributions  and  the  volume  frac¬ 
tions  of  YSZ  and  LSM,  an  increase  in  the  function  F3(I<)  leads  to  an 


Fig.  1.  Plot  of  F(a )  versus  a. 
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Fig.  2.  Plot  of  F2(0p)  versus  6>P. 


increase  in  LTPB.  Therefore,  it  is  of  interest  to  determine  if  there  are 
particle  shape(s)  that  increase  the  value  of  F3(I<)  substantially  above 
36.  For  such  an  analysis  it  is  imperative  to  use  the  same  geomet¬ 
ric  measure  of  “size”  for  particles  of  different  shapes,  because  the 
objective  is  to  determine  which  shape(s)  increase  the  LTpb  for  con¬ 
stant  values  of  mean  sizes,  size  distributions,  and  volume  fractions 
of  the  phases.  For  convex  particles,  the  orientation  averaged  par¬ 
ticle  caliper  diameter8  is  a  rigorously  defined  unique  measure  of 
“size”  that  can  be  used  to  compare  the  “sizes”  of  convex  particles  of 
different  shapes  (say,  plates  and  spheres).  For  a  spherical  particle, 
the  caliper  diameter  is  the  same  in  all  orientations,  and  therefore, 
the  orientation  averaged  caliper  diameter  of  a  sphere  is  equal  to 
its  diameter.  The  orientation  averaged  caliper  diameter  D  of  any 
convex  particle  can  be  computed  by  using  Minkowski’s  equation  of 
integral  geometry  given  below  [24]. 


(41) 


In  Eq.  (41 )  H  is  the  local  mean  curvature  of  an  infinitesimal  sur¬ 
face  element  dA  on  the  smooth  surface(s)  of  the  convex  particle  and 
the  surface  integral  is  to  be  performed  over  all  smooth  surfaces  of 
the  particle.  The  second  term  arises  from  the  edges  (if  present)  on 
the  particle  surface(s)  and  x  is  the  dihedral  angle  of  an  edge  ele¬ 
ment  of  length  dX.  The  second  integral  is  to  be  performed  over  all 
the  edges  of  the  particle.  If  the  particle  has  only  smooth  surfaces 
(for  example,  ellipsoids),  then  the  second  integral  is  zero.  In  the 
particle  population  having  a  distribution  of  D  values,  the  mean  size 
is  denoted  by  (D).  In  Eq.  (34),  (Dy)  and  (DL)  are  these  measures  of 
the  mean  sizes  of  the  YSZ  and  LSM  particles. 

For  parametric  analysis  it  is  convenient  to  use  cylinder  as  a 
model  shape  to  generate  equiaxed,  plate-like  (or  flake-like),  and 
needle-like  morphologies  by  varying  the  ratio  of  the  cylinder  radius 
R  and  length  I.  The  volume  Vcyl,  surface  area  Scyl  of  a  cylinder  are 
given  as  follows. 


Applying  Eq.  (41 )  to  surfaces  and  edges  of  a  cylinder  yields  the 
following  result  for  the  orientation  averaged  caliper  diameter  D  of 
a  cylinder. 


(44) 


Combining  Eqs.  (42)-(44)  yields  the  following  result. 


^cyl  = 


WL/R)  D3 

[JT  +  (L/R)f 


and, 


^cyl  — 


8zr[  1  +  {L/R)]q2 
[tt  +  (L/R)}2 


(45) 


(46) 


For  cylindrical  YSZ  and  LSM  particles  having  constant  (R/L)  ratio, 


one  can  write: 

,,  8jt{L/R)  o 

<Vy>  = - k-!—^(D1) 

(47) 

8„|,+(1/R)| 

[n  +  {L/R)\2 

(48) 

TO-  S’'a/R\<D=, 

[7t  +  (L/R)f 

(49) 

(Sl,.84.+(W1(d;) 

[7 r  +  (L/R)]2 

(50) 

(Dy),  (D3),  (D2),  (Dy)  are  the  population  average  values  of  D3  and 
D2  for  YSZ  and  LSM  size  distributions,  respectively.  Comparing  Eqs. 
(47)-(50)  with  Eqs.  (22)-(25)  yields  the  shape  factors  KjY,  Kil,  f<2Y> 
I< 2L-  Substituting  these  shape  factors  in  Eq.  (27)  gives  the  following 
result  for  cylindrical  particle  populations  of  YSZ  and  LSM  having 

constant  (R/L). 

F3(K)=[\+(L/R)]2[\+7r(R/L)]2 

(51) 

Vcyl  = 


Scyl  —  2tT 


(42) 

(43) 


8  Caliper  diameter  is  equal  to  the  distance  between  two  parallel  tangent  planes 
of  a  particle.  For  a  convex  particle,  there  are  two  (and  exactly  two)  parallel  tangent 
planes  for  every  angular  orientation.  Therefore,  there  is  a  unique  caliper  diameter 
for  each  orientation.  Caliper  diameter  can  vary  with  the  tangent  plane  orientation. 
The  orientation  averaged  caliper  diameter  D  is  obtained  by  averaging  the  caliper 
diameter  over  all  tangent  plane  orientations. 


Fig.  3  shows  a  plot  of  F3(I<)  versus  ( R/L ).  Observe  that  (i)  for  any 
value  of  (R/L),  F3(I<)  for  cylindrical  particles  is  higher  than  that  for 
spherical  particles,  (ii)  F3(I< )  reaches  high  values  for  large  ( R/L )  ratios 
that  correspond  to  plate-like  or  flake-like  particle  shapes,  and  (iii) 
F3(/C)  reaches  high  values  also  when  (R/L)  approaches  zero,  which 
corresponds  to  needle-like  shapes.  Even  for  equiaxed  cylindrical 
particles  ((R/L)  =  0.5),  F3(I<)  =  59.44,  which  is  higher  than  the  value 
of  36  for  the  spherical  particles.  For  (R/L)  equal  to  2,  F3(I<)  is  equal  to 
119.25:  an  increase  of  more  than  a  factor  of  three  over  the  value  of 
36  for  spheres.  Eq.  (34)  reveals  that  for  any  given  YSZ  and  LSM  size 
distributions  and  volume  fractions,  Ltpb  is  directly  proportional  to 
F3(I<).  Thus,  it  can  be  concluded  that  for  given  volume  fractions  of 
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Fig.  3.  Plot  of  F3(I< )  versus  K. 


the  phases  and  particle  size  distributions,  LTPB  depends  significantly 
on  the  particle  shape,  and  it  is  the  lowest  for  spherical  particles. 
Therefore,  changing  the  YSZ  and  LSM  particle  shapes  from  spherical 
to  plate-like,  or  flake-like,  or  needle-like,  can  substantially  increases 
Ltpb-  It  is  reported  in  the  literature  that  at  least  LSM  particles  can  be 
produced  in  a  flake-like  morphology  [48].  Thus,  the  analysis  brings 
an  interesting  opportunity  to  optimize  LTpb  via  a  suitable  choice  of 
particle  shapes. 

4.  Summary  and  conclusions 

An  analytical  expression  is  derived  for  total  triple  phase  bound¬ 
ary  length  per  unit  volume  LTpb  in  the  microstructure  of  composite 
LSM/YSZ  cathode.  It  is  assumed  that  the  microstructure  is  isotropic 
uniform-random  and  the  microstrucutral  volume  of  the  composite 
cathode  is  at  least  of  the  size  of  the  smallest  representative  volume 
element  of  the  global  microstructure,  and  there  is  no  significant 
coarsening  during  the  sintering  process  (which  is  reasonable  for 
LSM/YSZ  composite  cathodes  sintered  below  1200  °C).  The  model 
explicitly  brings  out  the  effects  of  particle  shapes,  mean  sizes,  the 
coefficient  of  variation  and  the  skewness  of  the  size  distributions, 
and  volume  fractions  of  YSZ,  LSM,  and  porosity  phases  on  LTPB.  The 
parametric  analysis  leads  to  the  following  conclusions. 

1.  Ltpb  is  the  lowest  when  YSZ  and  LSM  particles  have  spherical 
shape,  and  it  can  be  substantially  increased  by  using  flake-like, 
plate-like,  or  needle-like  YSZ  and  LSM  particles. 

2.  Mono-size  YSZ  and  LSM  particle  populations  lead  to  higher  ITPB 
as  compared  to  the  populations  having  a  distribution  of  sizes. 
Size  distributions  with  large  spread  (high  CV)  can  reduce  the 
Ltpb  by  almost  an  order  of  magnitude.  For  a  given  CV,  size  distri¬ 
butions  with  negative  skewness  lead  to  a  higher  value  of  LTPB. 

3.  For  given  relative  proportion  of  YSZ  and  LSM,  and  given  mean 
sizes,  CV,  and  skewness  of  YSZ  and  LSM  size  distributions,  the 
highest  value  ofLTPB  is  obtained  at  the  lowest  value  of  the  porosity 
volume  fraction  that  yields  sufficient  percolation  and  connectiv¬ 
ity  of  pores  for  gas  flow. 

4.  In  the  microstrucutral  regime  of  interest  for  SOFC  applications 
where  all  the  three  phases,  LSM,  YSZ,  and  porosity,  must  topo¬ 
logically  percolate,  the  relative  proportion  of  YSZ  and  LSM  does 
not  significantly  affect  Ltpb- 
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